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Abstract 


We present two- and three-dimensional simulations of the growth of quantum 
fluctuations of the scalar fields in supersymmetric hybrid inflation models. For 
a natural range of couplings, sub-horizon quantum fluctuations undergo rapid 
growth due to scalar field dynamics, resulting in the formation of quasi-stable 
non-topological solitons (inflaton condensate lumps) which dominate the post¬ 
inflation era. 


1 mattb@amtp.liv.ac.uk 2 j .mcdonald@lancaster.ac.uk 



1 Introduction 


Hybrid inflation models [T] are a promising class of model, having a classically flat 
inflaton potential without requiring fine-tuning of coupling constants. In order to 
control radiative corrections to the inflaton potential, supersymmetric (SUSY) hybrid 
inflation models are favoured Email- D-terrn hybrid inflation models are of par¬ 
ticular interest, since they naturally account for the absence of order H 2 corrections 
to the inflaton mass squared term coming from non-renormalisable interactions with 
non-zero F-terms. However, in spite of requiring fine-tuning to eliminate order H 2 
corrections, several F-term hybrid inflation models have also been proposed. In light 
of recent reconsideration of the issue of fine-tuning in SUSY models and the idea of 
’split supersymmetry’ j3], such models may still play an important role. 

An important issue is the process by which inflation ends and the Universe reheats. 
It is known that the end of inflation and reheating in hybrid inflation is dominated 
by the non-linear dynamics of the scalar fields responsible for inflation PPP- Sub¬ 
horizon quantum fluctuations are rapidly amplified by scalar field dynamics once the 
symmetry breaking transition ending inflation occurs, becoming effectively classical 
modes PEj. The energy density at the end of inflation rapidly becomes dominated by 
spatial fluctuations of the inflaton. For the case of D-terrn hybrid inflation, the rapid 
growth of the perturbations is expected to dominate the energy density for A > 0.1 g 
PJ. Similar dynamics may be observed in models of tachyonic inflation [TO] . 1 

What is the subsequent evolution of this inhomogeneous scalar field? There are 
two possible outcomes. The growth of the spatial perturbations of the inflaton could 
result in scalar field waves which scatter from each other, dispersing the energy of the 
initially highly inhomogeneous state pm- Alternatively, the growth of the spatial 
perturbations might continue, resulting in the formation of non-topological solitons 
corresponding to inflaton breathers 0 (also known as oscillons (I2J). This was first 
suggested in an alternative view of the end of hybrid inflation, inflaton condensate 
fragmentation |S|, which was based upon perturbing a homogeneous inflaton con- 

1 For related analyses of tachyonic preheating see GH 
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densate assumed to form immediately after the end of hybrid inflation. In the 
non-topological solitons were called inflaton condensate lumps. Their existence can 
be understood as the result of an effective attractive interaction between the inflaton 
scalars due to the symmetry breaking held |S]. Although the assumption of a homo¬ 
geneous inflaton condensate is violated by the rapid growth of quantum fluctuations 
within a few coherent oscillation periods, the possibility that spatial perturbations 
could form non-topological solitons remains. Such objects were subsequently observed 
in a numerical simulation of the phase transition in hybrid inflation H3|. Oscillons 
such as inflaton condensate lumps are known to be quasi-stable, having a lifetime 
much longer than the period of the inflaton coherent oscillations [[[2]. They decay by 
a slow radiation of scalar held waves. For the case of the effective theory associated 
with SUSY hybrid inflation models, the lifetime of the inflaton condensate lumps is 
10 3 — 10 4 times the oscillation period 

The possibility of non-topological soliton formation at the end of SUSY hybrid 
inflation is strengthened by the existence of stable Q-ball solutions of the scalar held 
equations of D- and F-term inflation H3- These are two-held Q-ball solutions, com¬ 
posed of a complex inflaton held and a real symmetry breaking held. They carry a 
conserved global charge associated with the inflaton. As a result, it is possible that 
SUSY hybrid inflation could end by initially forming neutral non-topological solitons 
which subsequently decay to pairs of stable Q-balls of opposite global charge. The 
formation of neutral condensate lumps which decay to Q-ball pairs has been observed 
numerically in the context of Afheck-Dine condensates along MSSM hat directions [S]- 
(See also T7J.) It occurs because the initially neutral condensate lumps are unstable 
with respect to growth of perturbations of the phase of the complex scalar held, result¬ 
ing in the formation of an energetically preferred state consisting of a stable Q-ball, 
anti-Q-ball pair. This is likely to be a generic feature of the evolution of quasi-stable 
neutral condensate lumps whenever a related stable Q-ball solution exists. 

Thus there exists the possibility that the post-inflation era of SUSY hybrid inflation 
models will be dominated by Q-balls containing most of the energy of the Universe. 
Whether this happens depends upon whether the inhomogeneous state of the Universe 
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immediately following inflation evolves into quasi-stable neutral condensate lumps or 
simply disperses as scalar field waves. In this paper we report the results of a numerical 
simulation of the growth and evolution of quantum fluctuations of the scalar fields 
in SUSY hybrid inflation models. The growth of quantum fluctuations into classical 
fluctuations is studied using the equivalent classical stochastic field method JTSj, which 
replaces the quantum field theory by an equivalent classical field theory with a classical 
probability distribution for the initial conditions. We focus on the case of D-term 
inflation. The results for F-term inflation are likely to be similar to those of D-term 
inflation for a particular choice of the D-term inflation model couplings. 

The paper is organised as follows. In Section 2 we review SUSY hybrid inflation 
models and the scalar field equations. In Section 3 we discuss the initial conditions 
and details of the numerical simulation. In Section 4 we present the results of 2-D and 
3-D simulations of the evolution of the scalar fields. In Section 5 we discuss a scaling 
property of the D-term inflation held equations and the dependence of our numerical 
results on the coupling constants. In Section 6 we compare our results with earlier 
numerical simulations. In Section 7 we present our conclusions. In the Appendix we 
review the equivalent classical stochastic held method. 

2 SUSY Hybrid Inflation Models 

The minimal D-term inflation model is described by the superpotential |2| 

W = AS4>+4>_ , (1) 

where S is the inhaton, $-)- are helds with charges ±1 with respect to a Fayet-Iliopoulos 
U(1)fi gauge symmetry, £ > 0 is the Fayet-Iliopoulos term (£ 1//2 = 8.5 x 10 15 GeV |EJ) 
and g is the U(1)fi gauge coupling. The scalar potential is given by 

V = A 2 1-S'! 2 (|$ + | 2 + |<F_| 2 ) + A 2 |<F + | 2 |<F_| 2 + 9 — (|<F+| 2 - |<f>_| 2 + £) 2 . (2) 

Once | S'| < s c / a/2 = g^ 1 ^ 2 / A, <f>_ develops a non-zero expectation value, breaking the 
U(1)fi symmetry and ending inflation. 
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The simplest F-term inflation model is described by the superpotential jB] 

W = |s («> 2 - y) , (3) 

where S is the inflaton and $ is the field which terminates inflation, fi 2 and ij may be 
chosen real and positive. The scalar potential is then 

V = r/ 2 |S| 2 |$| 2 + ^|<F 2 -/i 2 \ 2 . (4) 

For the case of real scalar fields, the F-term scalar field equations become equivalent 
to the D-terrn equations in the limit where $ + = 0 and A = \/2 g. 

The superpotentials Eq. Q and Eq. m are invariant with respect to an R-symmetry 
under which only S transforms, which manifests itself as a global U(l) symmetry of the 
scalar potential, S —> e ia S. It is this symmetry which is responsible for the conserved 
global charge of the hybrid inflation Q-balls HE). 

In the following we will focus on the minimal D-terrn inflation model. We expect the 
end of inflation in F-term inflation models to be similar to D-term inflation in the case 
A = y/2g, up to the effect of U(1)fi symmetry breaking and cosmic string formation, 
which our results indicate are of secondary importance to the evolution of the energy 
density. During inflation, the inflaton evolves due to the 1-loop effective potential, 
AV’dS'l, |$_|) 2 . This determines its rate of rolling at the end of inflation, an important 
initial condition for the subsequent non-linear growth of the field perturbations. We 
assume throughout that $ + = 0. This is justified since the $ + field has a large positive 
mass squared both during and after inflation. The equations of motion are then 

S + 3HS-^S=-\ 2 \$_\ 2 S-^ , (5) 

and 

+ 3= -AV| 2 $- + g 2 (f - |<F_| 2 ) . (6) 

2 AF(|5'|, |$_|), has been calculated in j5], where AV is given for s = V^Re(S) and <p- = \^2Re^-. 
The complex version is given by replacing s, b- by v / 2|S'|, v / 2| < t > _| in A V(s,<j>-). 
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3 Initial Conditions and Numerical Methods 


The mechanism behind the formation of spatial inhomogeneities of the inflaton field 
at the end of D-term hybrid inflation is the dynamical amplification of sub-horizon 
quantum fluctuations of the scalar fields. In the Appendix we review the formation of 
classical scalar field perturbations from dynamical amplification of quantum fluctua¬ 
tions. The process can be summarised as follows. At the end of hybrid inflation, the 
symmetry breaking transition results in the formation of an effectively tachyonic scalar 
potential for the inflaton. Rapid growth of the quantum modes in this potential results 
in semi-classical quantum fluctuations. If the quantum fluctuations grow to become 
semi-classical whilst the field equations are linear in the inflaton, then the quantum 
field theory becomes equivalent to solving the classical field equations with a classical 
probability distribution for the initial conditions, given by the Wigner function of the 
field modes and their conjugate momenta in the semi-classical limit. The subsequent 
evolution of the scalar field fluctuations can then be studied by solving the classical 
scalar field equations with stochastic initial conditions (equivalent classical stochastic 
field). In this way we can simulate the growth of sub-horizon quantum fluctuations at 
the end of hybrid inflation. 

3.1 Initial Conditions 

The initial conditions for the simulation are obtained from the values of 0(f) and 0(f) 
for the homogeneous zero modes and values of the equivalent classical momentum 
modes q( k, f) and p(k, f). (Here 0 represents any real scalar field.) The classical initial 
conditions are applied at the end of hybrid inflation, |Sj = s c /y/ 2, at which time the 
S and <f>_ fields are massless. 

The fields are studied in a comoving box of volume V = L 3 with periodic boundary 
conditions, with the classical field 0 and its derivative, 0, being expanded in terms of 
modes as 

^.( X ,«) = -)=j: 5 (k,«)e it - (7) 
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and 


0(x, t) = -^=J2 p( k > *) e<k ' x > ( 8 ) 

where k = Afc(a, / 5, 7 ) with integer a, /3, 7 and Ak = 27 t/L. 

For a massless scalar held in de Sitter space, the initial condition for non-zero 
modes q(k,t) at t — 0 corresponds to a complex Gaussian distribution for |g(k, 0)| 
with root mean squared (r.m.s.) value given by (Appendix) ( IS] 

|g(k,°)| r . ms = - 7 =( 1 + ^- i ) , (9) 


and a random phase factor e* 5k for each k. Here rj is the conformal time at the end 
of inflation (7 = —H^ 1 ) and k — |k|. Since </>(x, t) is a real scalar held, the modes 
and phase satisfy q(— k, t) = <p(k, t) and 5_k = —5k- The initial value of the mode 
function p(k, t) has an r.m.s. magnitude 

Wk,0)U * = 72%j( 1 + iv) ' (10) 


For the sub-horizon modes of interest, \kr]\ = \k/H\ 1. In this case, from Eq. Q and 

Eq. ([Top , initially |p(k, 0)| = |g(k, 0)| ~ H\q(k, 0)|. Since the mode g(k, t), correspond¬ 
ing to a wavefunction of wavenumber k, will satisfy g(k, t) ~ kq(k,t) H\q(k, 0)| at 
subsequent times, the initial value of q(k,t) from Eq. (FTTH) will have a negligible effect. 

In our simulations we modelled the initial mode function q( k, 0) by the r.m.s. value 
of |g(k, 0) | multiplied by a random complex phase for each k. To simulate the evolution 
of the scalar held we considered a spatial lattice with x = A x(i,j, k ) for integer i,j, k. 
The sum over non-zero modes in Eq. CD then gives the initial spatial perturbation on 
the lattice, 


, cos ( AkAxai + 9 a ) cos ( AkAx/3j + 9p) cos (AkAx'yk + 0 7 ) , 

( 11 ) 

where k a ^ = Ak[c? + /3 2 + 7 2 ) 1 / 2 and we have substituted 5k — 5 ► 9 a + 0@ + # 7 , where 
0 a> p^ are random phases. We then solved the scalar held equations numerically on the 
lattice using Eq. m and 50(x, 0 ) = 0 for each of the real scalar holds. 


<5<K X ,0) = 
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Only those modes which are amplified, such that their occupation number be¬ 
comes large, can be considered classical. Modes which are not amplified remain as 
quantum fluctuations. In simulations the quantum fluctuations must be regularised, 
since summing over large k modes will result in a large initial amplitude for the spa¬ 
tial perturbations, causing the assumption of linearity of the held equations to break 
down. We introduce a cut-off A at a value given by the mass of the inflaton held in 
vacuum, A = 777,5 = A^ 1 / 2 . This should be close to the largest value of k for which 
modes can grow dynamically, k max , although the exact value of k max cannot be known 
a priori since it is determined by the full scalar held dynamics. 

4 Results of Numerical Simulations 

In order to investigate the formation of neutral inflaton condensate lumps, we have 
studied a model with a real inflaton held, s = y/2Re(S). We considered numerical 
simulations in two- and three-dimensions, using a staggered leapfrog routine and the 
initial conditions for quantum fluctuations Eq. <HB applied to s, (ft- 1 and (ft _ 2 at 
s = s c , where <l>_ = {(ft- 1 + icft- 2 )/\/ 2 . 

In Figure 1, we show the results for the growth of the energy density in a 2-D 
simulation using a 200 x 200 lattice, for the case A = 0.14 and g — 1, for four different 
times. We compare the energy density p with the mean energy density averaged over 
the lattice, p. The time from the hrst to last hgure corresponds to hve S coherent oscil¬ 
lation periods, Ts, where t$ = 27T m^ 1 . The growth of the energy density perturbations 
to non-linearity occurs within the hrst two coherent oscillation periods. We find that 
roughly circular lumps of energy density form, which are stable over the duration of the 
simulation. Essentially all of the energy density becomes concentrated in the lumps. 
We also observe some condensate lumps coalescing into larger lumps, in particular by 
comparing Figures lc and Id. The lattice spatial dimension is given by L — 207rm^ 1 , 
which is small compared with the horizon H~ l « ( 3 / 47 r) 1 / 2 A 2 Mp;/(pm|). 

In Figure 2, we show the amplitude of the S held at the centre of one of the 
energy density lumps in Figure 1. The oscillations of an unperturbed homogeneous 
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(a) t = t s 


(b) t = 3 t s 




(c) t = 4 T S 


(d) t = 5 t s 


Figure 1: 2D growth of energy density perturbations for A = 0.14 and g = 1 on a 200 
x 200 lattice. Black areas correspond to p > 4p. 
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inflaton field are shown for comparison. The field inside the lump is oscillating in time, 
showing that the inflaton condensate lump is indeed a quasi-stable oscillon. This is 
important, as it is not obvious that the energy density lumps in Figure 1 are not 
simply concentrations of the energy density with no coherent structure. If this were 
true, it would be possible for such concentrations to have a short lifetime due to a 
large amount of kinetic energy associated with their constituent scalar particles. The 
fact the energy density lumps are oscillons suggests that they will have a very long 
lifetime compared with the coherent oscillation period, as observed in numerical and 
analytical studies of single-field oscillons mm- 



Figure 2: Coherent oscillations of an energy density lump. 

In Figure 3 we show the results of a 3-D simulation on a 50 3 lattice. The regions 
where the energy density is greater than 4 p are shown. This clearly demonstrates the 
formation of inflaton condensate lumps in the 3-D case, confirming the results of the 
higher resolution 2-D simulations. 

We find no evidence of cosmic strings forming within the volume of our our lattice. 
This is consistent with our previous semi-analytical analysis, which indicated that the 
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Figure 3: 50 3 lattice 3D simulation. The plot shows regions of energy p > Ap at t = 5 ts- 
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typical seperation of U(1)fi cosmic strings formed at the end of hybrid inflation would 
be not very much smaller than the horizon, corresponding to the wavelength of the 
first <f>_ fluctuations to be able to grow classically and break the U(1)fi symmetry [9]. 


5 A Scaling Property of the D-term Inflation Field 
Equations. 


Throughout we have considered the case g = 1. However, by rescaling the spacetime 
coordinates, the D-term inflation held equations can be written in a form that depends 
only upon the ratio X/g. As a result, the condition for the formation of inhaton 
condensate lumps is likely to depend primarily upon the ratio X/g, with only a weak 
logarithmic dependence on g due to the initial value of the perturbations. 

Under the rescaling of the coordinates t = t/g, x = x/g, the tree-level scalar held 
equations become 

V 2 A 2 

S + 3HS —-—S' = —-|<h_| 2 S, (12) 

a 2 g 2 


and 


V 2 A 2 

<h_ + 3 M_ - — = — -\s | 2 <h_ + (/ 

a 2 g 2 v 


$_|U <h_ , 


(13) 


where derivatives are now with respect to x and t. We have defined H = a~ l da/dt = 
H/g. H is a function purely of X/g. This follows since H 2 oc p where, in terms of 
derivatives with respect to rescaled spacetime coordinates, 


p = g 2 |A| 2 + |<i>_| 2 + 


VS 


+ 


Vd>_ 


+ V(S,$-)/g 2 \=g 2 p 


(14) 


V(S, $_)/^ 2 is purely a function of X/g, therefore p depends only on X/g. Thus 
H oc p 1 / 2 = gp 1 / 2 and so H — H/g oc p 1 ? 2 . Therefore H = H(X/g). 

The dependence of the rescaled held equations only on X/g does not extend to the 
quantum fluctuation initial conditions, Eq. (ED. In terms of rescaled coordiantes and 
the rescaled wavenumber, k = k/g, the initial fluctuation becomes 


1 


¥(x,°) = 


8 g 


- 1/2 


2k 


cos ( AkAxai + 0 cos ( AkAxfdj + 6p) cos ^ AkAx'yk + 6 , 


aft') 


(15) 
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The upper limit of the sum, corresponding to the largest value of k for which the 
modes can grow, is determined from the rescaled field equations and is therefore purely 
a function of X/g. Thus the initial conditions for the rescaled field equations are 
proportional to g -1 / 2 . However, the growth of perturbations in a tachyonic potential is 
typically rapid (exponential), in which case the condition for perturbations to become 
non-linear will depend only weakly (logarithmically) on the initial perturbation. As 
a result, if inflaton condensate lumps form for one value of g (we considered g — 1 
in our simulations) then they will form for any value of g for the same value of X/g, 
up to a weak dependence on the initial conditions. The effect of reducing g will be 
to increase the time for the lumps to form (t = t/g) and to increase the size of the 
resulting lumps {x = x/g). 

6 Comparison with Previous Simulations 

In |Ej a simulation of SUSY hybrid inflation was reported for the case equivalent to D- 
terrn inflation couplings A = 0.014 and g = 0.14, such that X/g = 0.1. It was concluded 
that the energy density would be transferred to inhomogeneous classical scalar field 
waves. Although no non-topological soliton formation was observed, this value of X/g 
is on the borderline of values for which we would expect to find non-topological soliton 
formation 0. Therefore the results of ;H| are not in any obvious way inconsistent with 
the results we have obtained here. 

In [23 a simulation of the end of D-term was performed, again for the case equiv¬ 
alent to A = 0.014, g = 0.14 in D-term inflation. The formation of non-topological 
solitons was observed along the length of cosmic strings. However, only quantum 
fluctuations of the symmetry-breaking field were considered in pi, with no inflaton 
fluctuations, which cannot be considered fully realistic. As discussed above, we expect 
that fluctuations of the symmetry-breaking field will only play a minor role in the 
evolution of the total energy density, being primarily associated with cosmic-string 
formation on scales not very much smaller than the horizon |9]. Nevertheless, the 
results in j!3| confirm the possibility of non-topological soliton formation at the end 
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of SUSY hybrid inflation. 


7 Conclusions and Discussion 

Our numerical simulations have shown that growth of quantum fluctuations of the 
scalar fields in SUSY hybrid inflation models results in the formation of inflaton con¬ 
densate lumps corresponding to spherical lumps of coherently oscillating scalar field, 
also known as oscillons. The energy density of the Universe becomes almost entirely 
concentrated in these condensate lumps. We have focused on the case of D-terrn hy¬ 
brid inflation, in which case the inflaton condensate lumps form if A ^ O.lg. Although 
a full analysis remain to be done, similar results may be expected in the case of F-term 
inflation, whose classical field equations are similar to those of D-terrn inflation for the 
case A = V2g. 

We have performed simulations for the case of a real inflaton field. The fact that the 
objects which form in this case are oscillons is significant, since from numerical and an¬ 
alytical studies of single-field oscillons such objects are known to be quasi-stable, with 
very long lifetimes compared to the scalar field coherent oscillation period. Therefore 
we expect that the post-inflation era will consist entirely of long-lived inflaton conden¬ 
sate lumps, even though our present numerical simulations are unable to resolve the 
evolution over more than about 10 coherent oscillation periods. 

In the case of a complex inflaton field, oscillons are expected to be a precursor to 
stable Q-ball formation via decay to Q-ball pairs. Decay to Q-ball pairs is likely to 
occur long before the neutral inflaton condensate lumps decay via radiation of scalar 
field waves. Thus if inflaton condensate lumps form in the case of a real inflaton field, it 
is very likely that decay to Q-ball pairs and a Q-ball dominated post-inflation era will 
occur once a complex inflaton is considered. Since the Q-balls are classically stable, this 
era will only end once the scalar particles forming the Q-ball decay perturbatively to 
MSSM particles, reheating the Universe. A highly inhomogeneous inflaton condensate 
lump or Q-ball dominated era would have significant consequences for post-inflation 
SUSY cosmology. For example, it has been shown that the dynamics of scalar fields 
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in the post-inflation era will not be subject to order H mass corrections ra, with 
important consequences for Affleck-Dine baryogenesis, SUSY curvatons and moduli 
dynamics. Reheating will also be inhomogeneous, via Q-ball decay. 

We plan to develop the capabilities of our numerical simulations in the future, with 
the ultimate goal of achieving a complete understanding of the post-inflation era of 
SUSY hybrid inflation models. 
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Appendix. Non-linear Growth of Quantum Fluctu¬ 
ations and the Equivalent Classical Stochastic Field 

In this Appendix we review the growth of quantum fluctuations in the presence of a 
time-dependent tachyonic mass squared term. We first consider a real scalar field in 
Minkowski space with a positive or negative mass squared term. We then generalise to 
the case of a Friedmann-Robertson-Walker background, in particular de Sitter space. 

In order to derive the classical initial conditions used in the lattice simulations, 
we quantize the theory in a box of side L. We consider an initial state corresponding 
to a Minkowski or Bunch-Davies vacuum state for a massive scalar field. The scalar 
field theory is exactly quantized in terms of creation and annihilation operators for 
momentum modes, which is valid so long as interaction terms can be neglected in 
the Hamiltonian, or equivalently the scalar field equations are linear. To follow the 
evolution of the quantum fluctuations into the non-linear regime, the quantum field 
theory is replaced by an equivalent classical field theory with a classical probability 
distribution for the initial conditions EMU, equivalent in the sense that the results 
of quantum field theory for the expectation values of products of fields are equal to the 
values obtained by solving the classical field equations with the classical probability 
distribution for the initial field values. This can be done if the occupation numbers of 
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the quantum modes become large compared with one during the linear regime. The 
classical probability distribution is then given by the Wigner function of the held modes 
in the semi-classical limit. Although the classical probability distribution must become 
valid at some time during the period when the theory can be considered linear, once 
it becomes established it can be applied at any earlier time during the linear regime. 
The Lagrangian and Hamiltonian for the scalar held in Minkowski space are 

1 M )- 

C = -d^d^cj) - — 0 2 , (A-l) 

and 

H = 7T0 — C = ^ 7 r 2 + (V0) 2 + m 2 0 2 , (A-2) 

where the canonical momentum is n = dC/dcj) = 0. We consider the scalar held 
in a cube of volume L 3 with periodic boundary conditions. The held and canonical 
momentum are expanded in terms of momentum modes as 

<K X > t) = -J= Y, q{ k, t)e zk * (A-3) 

and 

<x,t) = -A^ p (k, i )e‘ t -, (A-4) 

where the components of k satisfy ki = 2nUi/L with n t an integer. 7r(x, t) = 0(x, t ) im¬ 
plies that p(k, t) = q (k, t). The non-zero equal-time canonical commutation relations 
for the held operators are 

[7r(x, t), 0(x , t)] = — id 3 (x — x ^ . (A-5) 

The held and conjugate momentum mode operators must then satisfy 

[p(- k > t), q{ 1, t)} = -iS K i . (A-6) 

In terms of mode operators the Hamiltonian becomes 

H = i^p f (k,t)p(k,t) +u'lq\k,t)q(k,t) , w£ = k 2 + m 2 , (A-7) 

2 k 
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where the mode operators satisfy the classical Hamiltonian equations of motion. In 
general, a Hamiltonian which is quadratic in the mode operators may be exactly 
quantized by expanding the operators in terms of time-independent creation and an¬ 
nihilation operators 

<?(k, t) = / k a k + / k al k (A-8) 


and 


p( k, t) = -igkCLy, + ig k al k 


(A-9) 


where [a k , oj] = 5 k i. p(k, t) = g(k, t) implies that g k = i/ k . The mode functions / k 
and g k satisfy the classical held equations for the momentum modes (which are, by 
assumption, linear in / k ) and the canonical commutation relations. The latter implies 
that the mode functions satisfy the condition (? k / k +(? k / k = 1. For the case of a constant 
positive mass squared term, the mode functions / k must also satisfy the condition that 
the Hamiltonian is correctly quantized, such that the energy eigenstates correspond 
to numbers of quanta of energy cn k . The unique solution for / k which satisfies both of 
these conditions is 


/k = 


0 -iu k t 


(A-10) 


To study the evolution of the mode functions from an initial time t — 0 at which 
m 2 > 0, we use Eq. (lA-lOj) as the initial condition for the mode functions. The 
subsequent evolution may then be studied by solving the classical equation of motion 
for / k so long as we are in the linear regime, since in this case the held equations may 
be satisfied by held operators expanded in terms of time-independent creation and 
annihilation operators. 

As an example, we consider the evolution for a model where a held is massless at 
t < 0 and gains a constant negative mass squared term once t > 0. The held operator 


at t > 0 satishes the equation, for modes with k 2 < \m\ 2 , 


q (k, t) = u>lq (k, t) , u; k = \m\ 2 - k 2 , (A-ll) 


with general solution 


g(k, t) = a k e^ kt + d k e u)kt 


(A-12) 
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and 


Ace-“ kt ) , 


(A-13) 


p( k, t) = q( k, t) = tn k (cpcA 1 ^ - 


where ct k and Ac are time-independent operators corresponding to linear combinations 
of a k and al k . Matching these with the vacuum (m = 0) solution at t — 0 implies 


that 


a k 


i 

2y / 2cJ k 





(A-14) 


and 


where cn k 
with 


Ac — 


2 v / 2cn k . 


i , Wk \ 

1 + —— ) « k 

o; k 


icn k \ + 

- o_ 

w k / 


(A-15) 


k|. Therefore g(k, t) for t > 0 can be written in the form of Eq. (1A-8I) 


/k 


i 

2y / 2cn k 



— ) + 

/ 




(A-16) 


For u) k f S> 1, g(k, t) ~ and p(k, t) ~ cn k a k e“ kt , so that to a good approximation 

g(k, t ) and p(k, t) commute, corresponding to the limit of classical physics. 

In order to study the growth of quantum fluctuations beyond the linear regime, 
we need the classical stochastic initial conditions which reproduce the held theory 
expectation values in the semi-classical limit. The classical probability distribution is 
obtained from the Wigner function, 


W(q, p , t) = / < q ~l,t\ P \ q+ X -,t 


> 


(A-17) 


where p = pi + ip 2 , q = q i + iq 2 and x = X\ + ix 2 . The q and p represent complex 
generalised coordinates and conjugate momenta, which in our case correspond to g(k, t ) 
and p(k, t). p = \ ^ 0 (t) >< ij} 0 {t) \ is the density matrix of the Schrodinger picture state, 

|>, which evolves from the initial vacuum state. In general the Wigner function 
gives the probability distribution of the observable q at t after integrating over p 
and vice-versa. In the classical limit the Wigner function tends towards a classical 
probability distribution, f(q,p,t), which reproduces the quantum held theory results 
for expectation values of the held operators. Although f(q,p,t ) should be calculated 
in the semi-classical limit, it can be applied at any time during the linear evolution 
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of the field, in particular at earlier times, since it satisfies the classical equations of 
motion )2U]. In the case of hybrid inflation we can therefore apply the classical initial 
conditions at the end of inflation (s = s c ), corresponding to massless scalar fields. 

The annihilation operators can be expressed as 


a k 


2 Re(hgl) 



-rrP(M) 

ml 


(A-18) 


By definition a k |0 >= 0, where the Heisenberg picture initial vacuum state, |0 >, is 
equivalent to the Schrodinger picture vacuum state at t — 0. In order to obtain the 
Schrodinger wavefunction -0(g, t) =< q\il) 0 if) > corresponding to the state evolving 
from the initial vacuum state, we express the equation a k |0 > as a wave equation 
in the position representation p3). In terms of the time evolution operator U(t), 
a k |0 >= 0 can be expressed as 

zft 

rf[q( k, 0) + -4p(k, 0)]U\0 >= 0 , (A-19) 

9 k 


where q(k,t) = Wq(k,0)U. Since U |0 >= | i^ 0 (t) > is the time-dependent Schrodinger 
state vector and q( k, 0) = g k , p( k, 0) = p k are the Schrodinger representation gener¬ 
alised coordinate and momentum operators, with p k conjugate to g_ k = g k ([p k , g k ] = 
—i), the momentum operator in the position representation becomes p k = — id/dq k . 
The wave equation can then be written as 

(ok + ^o (®c, ql,t) = o . (A-20) 

V ,9 k d 9kJ 

This has the solution |T^ 

1>o{qk,q*,t) = N 0 exp ^-^y|g k | 2 ^ , (A-21) 

where N 0 = (v^l/kl) -1 - The normalized probability distribution is then 

p(?k ' &t) = ^T“ p (“S) ' (A_22) 


This gives the probability that the magnitude of the classical mode amplitude q( k, t) 
at t is |g k |. This is a complex Gaussian distribution, with r.m.s. value for |g k | equal 
t° |/k|- Then g(k, t) = \q(k,t)\e zSk , with h k a random phase. 
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The Wigner function is then given by 


W(q k ,p k ,t) = J j d '\' ) d ' , ] i ( l(PklXl+Pk2X2) ^o(Qk-x/2,t)^ 0 (q k +x/2,t) . (A-23) 


Integrating gives 


W{q k ,p k ,t) = -^exp - |/k | 2 Pk - rr 7 ? 9 k 


(A-24) 


where F k = Im(fffg k ). In the limit \f k \ 2 3 > 1, where the Wigner function tends 
towards the classical probability distribution, the Wigner function is vanishing except 


F k 

Pk \M 2qk 


(A-25) 


Since the classical probability distribution can be applied at any time during the linear 
regime, we can choose to apply the initial conditions at the end of hybrid inflation, 
such that the fields are massless. For the case of flat space, the mode functions with 


rn — 0 are given by 


? 9 k ~ ^ 


(A-26) 


where k — |k|. Thus F k = Im(f^g k ) = 0. Therefore the flat space classical initial 
condition for p(k, t) is p(k, 0 ) = 0 , whilst the flat space classical initial conditions for 
g(k, t) at t — 0 corresponds to a complex Gaussian distribution for \q(k, 0) |, with r.m.s. 


value 


|?(k, 0 ) \rms 


(A-27) 


and with a random phase factor e* 5k for each k. 


These results can easily be generalised to a scalar field in Friedman-Robertson- 
Walker spacetime. In this case the action is 

/ 1 777 ^ 

d 4 x^/=g -d^(f)d^(f) - —(f) 2 . (A-28) 

In terms of y — a(t)(f) and conformal time rj (dt = a[t)di 7 ), where a(t) is the scale 
factor, this becomes 


S = d A x -\y - y 


(Vy ) 2 m 2 a 2 2 


(A-29) 
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where ' denotes differentiation with respect to conformal time and d 4 x is now under¬ 
stood as an integral over comoving coordinates and conformal time. The conjugate 
momentum is 


dC / a 

* = w-t = y —y 

ay a 


(A-30) 


and the Hamiltonian is 


H 


d 3 x 


vr 2 (Vy ) 2 

2 + 2 


2 2 ' 

m a 9 a 

+ —^—y + —*y 

2 a 


(A-31) 


Expanding ?/(x, 7]) and 7 r(x, rf) in terms of modes as before (but now with y in place 
of t and y in place of 0) implies that 


H = -J2 
2 ^ 
L k 


p( k, v)*p( k , v) + ^k?( k > f7) + ?( k , v) + — (p( k , y)q\K v) + p f ( k > v)y(.K v)) 

CL 

(A-32) 


where = k 2 + m 2 a 2 and 


p(k,y) = q(k,rj) 


—q(Kv) 

a 


(A-33) 


This is quantized by expanding q and p in terms of creation and annihilation opera¬ 
tors according to Eq. (ED and Eq. (ED- The mode functions satisfy the classical 
equation of motion 


/k - -/k = - ( k2 + mV) /k , (A-34) 

a v 7 


with the requirement that the initial /k reduces to the Minkowski modes in the limit 
k^> H (Bunch-Davies vacuum). 

For the case of de Sitter space, H = constant , the mode equation becomes 




0 , 


(A-35) 


where ij = —1/aH. For a constant m 2 , the solution which tends to the Minkowski 
vacuum in the limit k H is 

AM = (-to,) , S = 9 - , (A-36) 


20 













where H /b is a spherical Hankel function of the first kind. 

The Wigner function analysis is as before except with g in place of t. / k is now 
given by Eq. (|A-36J) . whilst g k follows from p(k, g) = q (k, rj) — {a /a)q(k,g), which 
implies that g k = i(/ k — (a'/a)/ k ). The initial conditions we impose correspond to de 
Sitter space with a massless held, in which case u = 3/2 and so 


/k = I rr 1 v /3 h#3 1) (“M = 
2 


e ~ikr] 

V2k 


1 - 


kg 


Ik 


ft = \i if ikn . 


where we have used the expression for the Hankel function, 


= 


( 2 X 1 / 2 ■ 

i' 

(— ) e 

1 + - 

\7TXj 

X. 


(A-37) 

(A-38) 

(A-39) 


Thus F k = J/n(/ k g k ) = 1/(2 kg). The classical r.rn.s. values for the magnitudes of the 
mode functions are then 

\q(k } g)\ rms = |/ k | = --L (l + (A-40) 

and 

|p(k '’ ,)U ' = 73M( 1 + ^?) ' (A " 41) 

These may be used to model the initial conditions for the massless scalar fields at the 
end of hybrid inflation. 

In practice, we numerically solve the classical scalar held equations for 0(x, t ) rather 
than y{yi,g). In general, y(^g) = a0(x, t) and 7 r(x, g) = a 2 0(x, t). Choosing a = 1 
at the end of inflation (corresponding to choosing g = —H^ 1 at t = 0 ) implies that 
</>(x, 0) = y(x, g) and 0(x, 0) = 7 r(x, g) at the end of inflation. In this case the mode 
functions in Eq. (IA-40|) and Eq. (IA-41D will provide the initial conditions for the mode 
expansions of </>(x, t ) and </>(x, t ) in de Sitter space. 
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